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Effect of randomness in logistic maps 
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We study a random logistic map Xt+i = atXt[l — Xt] where at are bounded ( q\ < at < ( 72 ), random 
variables independently drawn from a distribution, xt does not show any regular behaviour in time. 

We find that xt shows fully ergodic behaviour when the maximum allowed value of at is 4. However 
(a?t—► 00), averaged over different realisations reaches a fixed point. For 1 < at < 4 the system shows 
nonchaotic behaviour and the Lyapunov exponent is strongly dependent on the asymmetry of the 
distribution from which at is drawn. Chaotic behaviour is seen to occur beyond a threshold value 
of qi ( 52 ) when 52 (qi) is varied. The most striking result is that the random map is chaotic even 

when q 2 is less than the threshold value 3.5699. at which chaos occurs in the non random map. 

We also employ a different method in which a different set of random variables are used for the 
evolution of two initially identical x values, here the chaotic regime exists for all q\ ^ 52 values. 


PACS numbers: 05.45.-a,05.90.+m,87.23,74.40.De 

I. INTRODUCTION 

Many natural phenomena show chaotic behaviour and 
a possible route to chaos is provided by nonlinear dynam¬ 
ical models. In the study of nonlinear dynamics, logistic 
map is a well known area of research [l|. Study of lo¬ 
gistic map is relevant in population dynamics ;2J, image 
encryption , electronic circuit, pseudo-random num¬ 
ber generators [6, Q etc. Mathematically, the logistic 
map is given by 

x t+ i = ax t (l - x t ). (1) 

The logistic map is perhaps the most simple and illustra¬ 
tive example of iterative maps of the form Xt+i = f{xt) 
which shows very complex dynamical behaviour. In this 
map, for 1 < a < 3 a nonzero steady value is reached. 
For larger a < a c = 3.5699. multiple fixed point val¬ 
ues with periodicity 2,4, 6,.etc. occur, while above a c , 

chaos occurs which is exhibited by a positive Lyapunov 
exponent [lj. 

In reality, dynamical phenomena e.g. population dy¬ 
namics, usually involve an amount of stochasticity. In or¬ 
der to incorporate this feature, we analyse the behaviour 
of the logistic map where a does not have a fixed value 
but is a random variable at., drawn independently from a 
distribution. In some earlier works |8H10}, such randomi¬ 
sation of a has been considered where the main issue was 
to check whether Xt —> 0 as t —> 00 such that the distri¬ 
bution of x at large times is a delta function at x = 0. It 
was concluded that under certain condition this is true. 

We use different distributions for choosing the control 
parameter a*. We choose a* to be bounded, i.e., qi < 
at < <72 with qi > 1 to avoid the fixed point x = 0 as 
obtained in the earlier cases SHU- 

We study the behaviour of the quantity A t defined as 

A t = (\x t — x' t \), (2) 

where Xt. and x' t denote the two different evolutions. A t is 
the average value obtained from different configuration as 
denoted by the angular brackets on the right hand side of 


ED- The evolution has been implemented in two different 
ways: in the traditional method (TM), one starts with 
two initially close values Xg and x' 0 to calculate A t as 
these two initial x values are evolved using the same set 
of a*. In the other method, different sets of at are used 
for initially identical values of x. The latter is called the 
“Nature vs. nurture” method following El- 

As usual, if A t grows (or saturates at a finite value 
as it cannot increase indefinitely) we conclude that the 
chaotic regime is reached. Varying qi and q 2 we identify 
such regions in both the methods. Apart from identifying 
the chaotic region, we are also interested in comparing 
the TM and NVN methods which have led to different 
results in interacting dynamical systems 0 Hi- 

One of the main objectives is to study whether univer¬ 
sal behaviour is observed when we vary the distribution. 
For this purpose, we have used a uniform distribution 
and a triangular distribution which has a peak at a value 
a p (?i < a p < < 72 ), from which the a values are chosen. 
For the triangular distribution a p = (<71 + < 72 )/2 signifies a 
symmetric triangular distribution. Deviation of a p from 
this value quantifies the asymmetry. 

In the next two sections we discuss the results and in 
section IV summary and discussions are presented. 

II. ERGODICITY AND CONVERGENCE 

We first report some general results which are in stark 
contrast to the non-random map. We use either a uni¬ 
form distribution or a symmetric triangular distribution 
for these studies. 

Let the initial value of x be chosen randomly from 
[0:1]. In the non-random logistic map, it is known that x 
has finite number of attractors in the nonchaotic regime. 
In the random case however x has an ergodic type be¬ 
haviour in the sense that it does not reach a fixed at¬ 
tractor and can assume values between x m in and x ma x', 
A x = ( x max - x min ) is nonzero for q x ^ q 2 . When q 1 
is fixed at 1, A x shows an increasing behaviour with q 2 . 
However when (72 = 4, A^, cs 1 for any value of q\ > 1 
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showing that the system becomes fully ergodic (fig. [2]). 
The fate of two independent evolutions depends on qi 
and q 2 ] however we note that the average value (x*) of 
all such evolution shows convergence and is non ergodic. 
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FIG. 2. (Color online) TM results: The variations of A x , Xmin 
and Xmax with q 2 when q\ = 1.0 (lower panel). The variations 
of A x , Xmin and Xmax with q\ when <72 = 4.0 (lower panel). 
Both the figures are for uniform distribution. 
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FIG. 1. (Color online) TM results: Two different evolutions 
of x(t) in the traditional method. Left panel: uniform distri¬ 
bution; Right panel: symmetric triangular distribution. The 
steady state values averaged over many configurations also 
shown (see also figs. Up Hr). 

It is known that for a fixed value of a, the non-random 
logistic map has a fixed point at x = which is stable 
below a = 3. For the random case, one can approximate 
theoretically, 

<*) = ^ ( 3 ) 

a 

where a denotes the average value. If we allow the dis¬ 
tribution of a to vary between qi and q 2 , we find that 
the value of (xt^oo) obtained numerically differs from 
( x } obtained using eq. [3] (Table QJ. We have shown that 
while the non ergodicity behaviour is present for both the 
uniform and symmetric triangular distribution, the devi¬ 
ation of (x) is larger for the uniform distribution com¬ 
pared to that in the symmetric triangular distribution 
from (xt~y oo). Also, ((x) — (xt_ >00 ))/(x) increases with 
92 — qi in the uniform distribution. However, this in¬ 
crease appears to be weaker in case of the symmetric 
triangular distribution. 

III. RESULTS FOR TM METHOD 

1. Non-chaotic regime 

We choose a pair of initial values of x which are slightly 
different and allow them to evolve as a function of time 
(eq. [2J . Typical evolutions for qi = 1 and <72 = 4 of 
two initially close x values are shown in figs, [lji, [ljo. Al¬ 
though x does not attain a fixed point value, we find 
that At indeed goes to zero in an exponential manner for 
chosen values of q\ and q 2 (e.g. for q\ = 1,172 =4; see 
figs. |3ji, EJd) signifying regular or nonchaotic behaviour. 
However, here we find that At varies in a nonlinear man¬ 
ner with Aq; precisely, At ex -\/Ao exp(Af), whereas for 


non-random logistic map it is A t oc Ao exp(At). The Lya¬ 
punov exponent A depends strongly on < 71 , <72 ; it shows an 
increase as q 2 — qi is decreased. The value of the Lya¬ 
punov exponent is A ^ —0.186 for uniform distribution 
for 1 < at < 4. 

At has been studied for asymmetric triangular distri¬ 
bution for 1 < at < 4 (Fig. [3jc) - Here the Lyapunov 
exponent remains constant upto a p ss 1.7 and decreases 
with a p for higher values (Fig. [3]:1) . This signifies that as 
a p increases, At vanishes in a slower manner. 


TABLE I. Table for steady state value of x. 


Distribution 

Theor. value 
(x) = 

Actual value 

>00) 

Uniform, 1 < a < 4 

0.600 

0.499 

Uniform, 3 < a < 4 

0.714 

0.636 

Symmetric triangular, 1 < a < 4 

0.600 

0.549 

Symmetric triangular, 3 < a < 4 

0.714 

0.650 


2. Chaotic regime 

When q 2 is a variable and q\ = 1, A t —> 0 as t —> 00 
for any value of < 72 - However increasing qi, we note that 
At^oo may reach a nonzero value, e.g, when q\ = 3 and 
q 2 = 4 signifying a chaotic behaviour. Typical evolutions 
for qi = 3 and q 2 = 4 of two initially close x values are 
shown in figs. 0^, |Dr. In this case, the damage saturates 
to a nonzero value (Figs. 0}:, |4ji main plot). We have 
used here either a uniform distribution or a symmetric 
triangular distribution for at- Saturation value of the 
damage A sat = A t ^oo has been studied for both the 
distributions for different values of q\ and < 72 - One can 
keep 172 variable and fix qi and observe the onset of chaos 
at a threshold value of <72 (Fig. [5] upper panel). Calling 
this threshold value q%, we note that q% is a function 
































3 



FIG. 3. (Color online) TM results: (a) Scaled value of dam¬ 
age At against time for uniform distribution, (b) scaled value 
of damage At against time for symmetric triangular distribu¬ 
tion, (c) At against time for the symmetric triangular distri¬ 
bution on [1,4] with different values of a p and (d) the variation 
of Lyapunov exponent with the value of a p for the asymmetric 
triangular distribution on [1,4] peaked at a p . 


FIG. 4. (Color online) TM results: (a) and (b) Two differ¬ 
ent evolutions of x(t) in the traditional method in the chaotic 
regime. The steady state values averaged over many config¬ 
urations also shown, (c) and (d) At against time for qi = 4. 
Main plot for qi = 3 shows saturation of A* at a nonzero 
value. Inset for q\ = 2 shows data when A sa t —» 0. Left 
panel: uniform distribution; Right panel: symmetric triangu¬ 
lar distribution. 


of qi and decreases with q\ which is expected. This is 
true for both distributions (Fig. [5] upper panel inset). 
We note that the minimum value of q\ for the onset of 
chaos is ~ 2.6 for the uniform distribution and ~ 2.7 for 
the symmetric triangular distribution. Even for q -2 < a c , 
where a c is the threshold value for onset of chaos in the 
non random map, one can observe a chaotic region for 
both the distributions. Similarly one can keep q\ variable 
and fix <72 and observe the onset of chaos at a threshold 
value of qi (Fig. [5] lower panel). Calling this threshold 
value qf, we again note that it is a function of <72 (Fig. [5] 
lower panel inset). What is striking is the presence of a 
peak at around 3.6 which is very close to a c . It is found 
that for <72 ~ 3.1, the minimum value of q\ required for 
chaos is ~ 3.0. Note that this is the value above which 
bifurcations start occurring in the nonrandom map. In 
fact for 3.1 < (72 < 3.6, the threshold value q{ is weakly 
dependent on (72 and remains ~ 3.0 in the entire region. 
Above q -2 = 3.6, however, smaller values of qi allow chaos. 


IV. RESULTS FOR THE NVN METHOD 

We next discuss the results for the NVN method. In 
this case, two copies of Xt, initially identical, evolve in¬ 
dependently, i.e., using different random numbers drawn 


from identical distributions. Here also the control pa¬ 
rameter is chosen from uniform and symmetric triangu¬ 
lar distributions. We find that independent of the values 
of qi and <72, the time evolved values of the two copies 
never converge as long as q\ 7^ <72. Both copies evolve 
with different values at all times (Fig. [H]); however, the 
configuration average of course is the same as that in the 
TM case. 

The damage A t as a function of time initially increases 
and then takes a steady value. The saturation value of 
A t is nonzero for all the different regions of control pa¬ 
rameter (Fig. 0 as is expected from the evolution of Xt- 
This is again true for both the distributions. 

The saturation value of the damage A sat has been 
studied. When <71 is fixed and the upper limit (72 = <7 
is varied (qi < q < 4), the value of A sat increases with 
q right from q > qi as there is no threshold value of the 
chaos (Figs. 0,0). 

If we keep q± = q as a variable (g > 1) and (72 fixed, 
the saturation value A sat shows an interesting behaviour. 
Up to <72 ~ 3.0, it decreases with q. At a critical value 
of (72 ~ 3.5, it shows a non-monotonic behaviour, with a 
sharp rise close to q = 3.0 before decreasing to zero at q = 
q 2 . Beyond this critical value of < 72 , the initial decrease in 
A sat becomes less prominent and it increases right up to 
q = <72 indicating there is a sharp discontinuity of A sa t 
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FIG. 6. (Color online) NVN results: Two different evolutions 
of x(t) in the method. Left panels: uniform distribution; 
Right panels: symmetric triangular distribution. The steady 
state values averaged over many configurations also shown. 


FIG. 5. (Color online) TM results: Saturation value of dam¬ 
age A sa t shows onset of chaos as q\ is fixed and q 2 varied 
(upper panel); q 2 is fixed and q\ varied (lower panel). Insets 
show variation of threshold values qj against q\ (upper panel 
inset) and gf against q 2 (lower panel inset). 


at q = q2 (Figs. [Ht, lEtO- This behaviour is like that in 
TM (Fig. [5] main plot), however the difference is, one has 
an onset of chaos in TM such that A sat increases from 
zero while in NVN, it increases from a nearly constant 
non-zero value. 

One can estimate the maximum possible value of A sat 
assuming two completely uncorrelated maps as: 




FIG. 7. (Color online) NVN results: Damage as a function 
of time for different regions of control parameter. Left panel: 
uniform distribution; Right panel: symmetric triangular dis¬ 
tribution. 


A est = 


/ / {x\ — X 2 ) 2 P{xi)P(x 2 ) dxi dx 2 

Jo Jo 


1 /*! 


( x\ — x\ — 2 x\X 2 ) dx\ dx 2 


0 -'o 


where P(x 1 ), P(x 2 ) denote distribution of x\ and X 2 - 
Assuming P{x\) and P{x 2 ) to be uniform, Ag St = 1/6. 
Therefore the expected value of A est = yT/6 « 0.41 for 
uncorrelated maps. A sa t is indeed less than A es t for both 
NVN and TM. 


V. SUMMARY AND DISCUSSION 

In summary, we have studied the behaviour of random 
logistic maps where the parameter a in eq. Q] is a ran¬ 
dom variable. Xt shows semi or fully ergodic behaviour 
for such maps, however (xt->oo) attains saturation val¬ 
ues which differ from the theoretical mean field values 
as given by eq. [3] The deviations are less in case of a 
symmetric triangular distribution as it has less variance 
(Table I). 

It is known that randomness in linear systems may 
give rise to chaos fl3j . In the present model, we can 
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FIG. 8 . (Color online) NVN results: (a) and (b) Saturation 
value of damage A sa t when q\ is fixed and q 2 varied (upper 
panel), (c) and (d) Saturation value of damage A aa t when 
52 is fixed and q\ varied (lower panel). Left panels: uniform 
distribution; Right panels: symmetric triangular distribution. 


identify nontrivial nonchaotic behaviour even with both 
randomness and nonlinearity. In the nonchaotic regime 
we find the unconventional behaviour A t ~ y/A q exp(At). 
Here, one can estimate the Lyapunov exponent A. We 
observe that A shows nonuniversality in the sense it shows 


strong dependence on the asymmetry of the distribution 
which may be quantified by \a p — ( qi + qr 2 ) /21 . 

Onset of chaos is noted in the traditional method at 
threshold values of <71 ( 92 ) which are dependent on q -2 
(qi). Minimum values for onset of chaos is found to be 
qi ~ 2.6 (for the uniform distribution) and q\ ~ 2.7 (for 
symmetric triangular distribution) when <72 = 4.0. On 
the other hand, for 52 , the corresponding minimum value 
is ~ 3.1 with qi ~ 3.0 for both the distributions. The 
most striking result is even when 52 < o. c , we obtain 
chaotic region for qi > 3.0 for both the distributions 

(Fig.®. 

In NVN, no threshold values of < 71 , (72 are obtained. 
Chaos occurs for all <72 > <7i- However A sat shows in¬ 
teresting variation with qi and <72 (Fig. 0- In general, 
A sat increases as 52 increases but A sat is not trivially 
dependent on (72 — q\ and hence we conclude that A sat is 
a nontrivial function of both q\ and ( 72 . 

Although the maps are random, certain effects of the 
nonrandom maps seem to be present in the TM results, 
e.g., the peak of q\ occurs at q 2 ~ a c and minimum value 
of q\ for (72 < a c is ~ 3.0 where bifurcation starts occur¬ 
ring in the non random case. Another point that needs to 
be mentioned is that it has not been possible to estimate 
Lyapunov exponents in the chaotic regime as saturation 
values are attained within very short times. 

As had been observed earlier [ll],[I2 , the TM and NVN 
methods yield completely different results. As in the case 
of damage spreading in opinion dynamics model [l 2 | , here 
too we find that the chaotic regime is obtained for any 
nonzero value of <72 — < 7 i in the NVN method. 
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